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We provide a detailed Dynamic Renormalization Group study for a class of stochastic equations 
that describe non-conserved interface growth mediated by non-local interactions. We consider ex- 
' plicitly both the morphologically stable case, and the less studied case in which pattern formation 

^SJ , occurs, for which flat surfaces are linearly unstable to periodic perturbations. We show that the 

latter leads to non-trivial scaling behavior in an appropriate parameter range when combined with 
O ' the Kardar-Parisi-Zhang (KPZ) non-linearity, that nevertheless does not correspond to the KPZ 

I universality class. This novel asymptotic behavior is characterized by two scaling laws that fix the 

critical exponents to dimension-independent values, that agree with previous reports from rmmerical 
simulations and experimental systems. We show that the precise form of the linear stabilizing terms 
does not modify the hydrodynamic behavior of these equations. One of the scaling laws, usually 
' (— I ' , associated with Galilean invariance, is shown to derive from a vertex cancellation that occurs (at 

^ ' least to one loop order) for any choice of linear terms in the equation of motion and is indepen- 

' dent on the morphological stability of the surface, hence generalizing this well-known property of 

, the KPZ equation. Moreover, the argument carries over to other systems like the Lai-Das Sarma- 

■ Villain equation, in which vertex cancellation is known not to imply an associated symmetry of the 

, equation. 
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I. INTRODUCTION 



' Many surface growth systems can be described by stochastic differential equations that capture the evolution of 
O ! the surface height h{r,t) at time t and point r on a d-dimensional substrate [H, In many cases these equations are 
^ O ^ ' derived from first principles 0, 0] , thus linking the small scale physics of the growth process with a coarse grained 
description of the surface height. Sometimes, this procedure ends up phrasing a physical problem in which interactions 
are short-range and local, into a stochastic equation in which effective interactions turn out to be long-range and/or 
non-local due to kinetic constraints. An example is provided by diffusion limited growth Q (or erosion Q) of 
an aggregate, in which the local growth velocity depends on the global morphology of the evolving cluster, due to 
' shadowing effects produced by the most prominent surface peaks over less exposed areas. Thus, a local transport 
mechanism such as diffusion induces long-range correlations among positions at the aggregate surface. Dynamically, 
this may result into a morphological instability in the case of a growing surface (protuberances grow faster, amplifying 
deviations from a flat interface), or morphological stability in the case of an eroding one (protuberances decay faster 
so that deviations from a flat interface are attenuated). 

In principle, the non-local morphologically unstable case is more interesting than the non-local morphologically 
stable case, since the latter is known to lead to trivial scaling behavior, in the sense that the corresponding scaling 
exponents can be obtained from simple dimensional analysis as applied to the corresponding linearized equation 7, 8] . 
On the other hand, morphological instabilities are key to pattern forming systems 9], non-local linear instabilities 
}^ ' having been the subject of extensive studies in the past. Celebrated examples include the MuUins-Sekerka (MS) 
I instability in solidificatio n fiol , t he SafFman- Taylor (ST) hydrodynamic instability in fluid flow or the Daerrius- 

Landau (DL) instability |13l4l6j occurring in the propagation of a premixed laminar flame, to cite a few. 

The stochastic equations that describe morphologically unstable surfaces typically combine linear terms that account 
for the selection of a typical wavelength associated with that one among the height field modes, h\^{t) = J^[/i(r, i)]k, 
that has the largest growth rate [oj , with additional non-linear terms that stabilize the exponential growth associated 
with the morphological instability, and prevent blow up the surface height field. Here, J- stands for space Fourier 
transform and k is wave vector. In the context of surface growth, the Kardar-Parisi-Zhang (KPZ) non-linearity 
(V/i)^ is expected whenever the interface evolves in absence of conservation laws 0, Thus, a natural example for 
a model of unstable surface growth subject to external fluctuations is given e.g. by 

dth = vV^h-]CV^h+^{Vhf +r], (1) 
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where rj is an uncorrelated Gaussian noise with zero mean and constant variance Ho, i.e., 

(,7(x, t)rKx', t')) = 2no 5(x - x') S{t ~ t'). (2) 

In equation ([IJ, j/ < 0, /C > 0, and A are parameters. A negative value for v implements the morphological instability, 
that is countered both by the biharmonic and by the nonlinear terms. Equation ([l} is the noisy Kuramoto-Sivashinsky 
(nKS) equation whose deterministic limit is a paradigmatic model for chaotic spatially extended systems, arising 
in a variety of physical contexts, like thin solid films, interfaces between viscous fluids, waves in plasmas and chemical 
reactions, reaction-diffusion systems, and others 9]. Actually, the large scale behavior of the nKS systems is in 
the KPZ universality class both in one and in two spatial dimensions |19l423j . This is a remarkable manifestation of 
renormalization, by which the "surface tension" parameter v flows from a bare negative value to a positive renormalized 
value at large scales, as initially conjectured by Yakhot thirty years ago [23 |. 
Indeed, the KPZ (or noisy Burgers [1^) equation reads 

at/i = j/V2;i+^(V/i)2+r/, (3) 

where > 0. The KPZ nonlinear term is invariant under translations in the growth direction (/i — > h + c, with 
c an arbitrary constant), in the substrate direction (x — > a; + a) and, besides, it is symmetric under rotations and 
inversion about the growth direction 1]. Its ubiquity has led to a huge amount of work in the past, the KPZ equation 
having become one of the central problems in non-equilibrium interface dynamics [l|, 0, [26| . Actually, based on 
Dynamical Renormalization Group (DRG) arguments KPZ scaling is generically expected at the asymptotic 

state for non-conserved growth systems. However, such a theoretical expectation has been too scarcely confirmed by 
experiments (28| . One possible explanation is that short-time morphological instabilities (of the type occurring e.g. 
in the noisy Kuramoto-Sivashinsky system) may hinder the observation of the asymptotic behavior [5|. For instance, 
for appropriate choices of parameter values in the nKS equation, KPZ scaling may become unobservable even in 
numerical simulations Moreover, some of the hypothesis made in the derivation of the KPZ equation may not 
be realistic for many experimental situations. In particular, complex processes may take place in the system, such as 
diffusion, advection in a fluid, elastic stress in the aggregate, etc. that break down the assumption of local interactions. 

Following the above ideas, a number of non-local versions of the KPZ equation have been proposed — see e.g. (29l - [3l| 
and further references in [28| — in order to account for scaling exponent values seen in experimental conditions that 
differ from the KPZ values. However, some of these are ad-hoc generalizations that lack a clear physical justifica- 
tion. Moreover, most of them consider the morphologically stable condition [e.g, v > in Eqs. ([ij, ([3])] which, for 
uncorrelated time-dependent noise, leads to scaling behavior associated with linear equations as mentioned above. 



A. A family of non-local interface equations 

There are other ways to generalize interface equations, that do have a consistent physical interpretation. For 
instance, one can generate additional terms on the rhs of, say, Eqs. ([T]), ([3]), if one enforces full compatibility of the 
equation with the symmetries of the system [l^. Otherwise, starting from a first principles formulation of a specific 
growth system in terms of, e.g., a moving boundary problem, one can derive explicitly the evolution equation for h{r, t) 
through projection of the whole bulk dynamics onto the interface, under suitable approximations. For non-conserved 
diffusion-limited growth, this has been accomplished for instance in (gI. Is^. For epitaxial growth systems, see a review 
in, e.g., [3]. Working along these lines, we have recently identified a whole family of non-local equations that is 



formulated as a single stochastic differential equation |33j . This equation has many relevant, well-known physical 



systems as particular cases. In Fourier space, the equation takes the simple form 

dthuit) = {~uk^ ~ /Cfc" - Uk^)h^{t) + ^^[(v;i)2]k + 77k(t), (4) 

where k denotes the magnitude of the d dimensional wave- number k and ?7k(0 ~ -^(^) with 77 as in equation ([2]). 
Here, /i, m, n, /C, and M are positive parameters such that < /x < 2, m > 2, and n > m, see below. We will consider 
both the morphologically unstable case in which v < 0, as in e.g. the noisy Kuramoto-Sivashinsky system, and 
the morphologically stable case v > 0. The term with coefficient Af represents higher order stabilizing mechanisms 
that are required to make sense of parameter choices in which fluctuations do not dissipate at small scales, i.e. 
—i/k'^ — /C/c™ — )■ -l-oo for k — >■ -l-oo, as in the /i — > to = 2 nKS-type limit. Still, one of the results of the present work is 
to explicitly prove the irrelevance of the term with coefficient Af for the asymptotic scaling. Note that, whenever fi is 
not an even integer, the corresponding linear term in equation Q is a fractional Laplacian acting on the height field, 

(-VY/^Mr) = c.,,Pv/^^j^^I^dr', 0<M<2, (5) 
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where PV denotes Cauchy principal value and Cd.i_i are appropriate numerical constants [3J, |35| . The Fourier repre- 
sentation of ([5]) is given by 

J-[(-V2)^/2/i] = fc^/ik- (6) 

Indeed, only if fi is an even integer can ^ be inverse Fourier transformed to get a local operator in real space, as 
in the noisy Kuramoto-Sivashinsky case. Otherwise, the fractional Laplacian corresponds to convolution of height 
differences with an algebraically decaying kernel. In the context of equilibrium systems with long-range interactions, 
equation (O corresponds to the so-called "weak" form of the latter, see (S^j for a review. For the morphologically 
stable case, relaxation operators with the same shape as ([5|) have also been advocated to account for long-range 
interactions in processes defined both on Euclidean and on fractal lattices, see e.g. 0, Isil, [s^l and references therein. 

For the morphologically unstable case, the linear dispersion relation that appears in equation (|4|), tJk — — i^fc''— /C/c™, 
is common to many celebrated pattern forming systems. Thus, for fi = 1 and m = 3 we get the MuUins-Sekerka 
or Saffman- Taylor dispersion relation, while the Darrieus-Landau and Kuramoto-Sivashinsky cases are obtained for 
/X = 1, m = 2, or for /i = 2, m = 4, respectively. Thus, Q allows us to study a large class of models in which a 
non-local morphological instability is coupled with the KPZ non-linearity, which includes a number of equations that 
are important on their own rights, but also local models as nKS equation that have been exhaustively studied in the 
literature. 

Apart from the noisy Kuramoto-Sivashinsky equation, we have the stochastic generalization of the Michelson- 
Sivashinsky (sMS) equation [s^ |3^ for /x = 1, m = 2, that provides a standard model to describe propagation of 
premixed laminar flames [l6| . Moreover, the sMS equation is the small slope approximation of a nonlinear equation 
found for growing interfaces controlled by ballistic transport [40], and has been also derived from first principles 
for reactive infiltration in porous media [41]. Likewise, an interface equation that we term MSKPZ, in which the 
Mullins-Sekerka dispersion relation (p — I, m = S) is combined with the KPZ nonlinearity, describes thin film growth 
by chemical vapor deposition (CVD) or by electrochemical deposition (in galvanostatic conditions) Q when the 
attachment kinetics of aggregating particles is fast compared to mean velocity of the growing interface, and has been 
recently employed to quantitatively describe CVD experiments [42]. Moreover, the linear dispersion relation appearing 
in has been found to provide a phenomenological description for experimental systems in which transport is of a 
diffusive nature, using fi as a fitting parameter [43| . 

In [33| . a numerical study has been reported for the family of equations in one and two dimensions, together 
with a preliminary DRG study for fixed to = 2 as a representative case. Particular attention was paid to the 
morphologically unstable case v < for which, strikingly, non-trivial (non-linear) scaling occurs. Thus, under such 
conditions simulations suggest that, for each dimension d, the universality class changes at ^ — ZKPz(rf)- Unstable 
equations characterized by /i such that 2:kpz('^) fJ- ^ '2 are all in the KPZ universality class (as for the noisy 
Kuramoto-Sivashinsky equation), while the scaling behavior for small values < ^ < ZKPzid) is characterized by 
dimension independent exponents that are different from those predicted by dimensional analysis. Here, ZKPzid) 
denotes the dynamic exponent of the KPZ class for the given value of d. For the definition of the critical exponents 
in these systems see section HIl 

In the present work, we perform a detailed DRG study of @ for arbitrary substrate dimensions and considering 
both the morphologically stable and unstable conditions. Although the perturbative DRG has well known limitations 
[25I [26| , it has recently shown a remarkable explanatory power for various types of systems that include conserved or 
non-conserved versions of the KPZ nonlinearity, under morphologically stable (see (4^ . [45 ^ and references therein) or 
unstable [i^ conditions. We expect relevant information to be also obtained from the DRG flow on the asymptotic 
properties of equation (|4|). Indeed, our present work extends previous analytical results in a number of important 
aspects: (i) By generalizing the linear dispersion relation and the nonlinearity appearing in the stochastic equation 
to more general forms, we prove that a vertex cancellation occurs to one loop order in our perturbative expansion, 
inducing a Galilean scaling relation at the non-linear fixed points found in the numerical simulations of [33]. The 
ensuing equation generalizes both the KPZ and the Lai-Das Sarma- Villain [13, 13 equations. Thus, explicit symmetry 
under a Galilean transformation, that occurs for the former equation but not for the latter, is seen not to be a necessary 
condition for the mentioned scaling relation to apply. Moreover, this scaling relation also holds at the asymptotic state 
numerically characterized for sufficiently small fi values in the morphologically unstable case of equation Q. (ii) We 
consider in detail the RG fiow for arbitrary d, obtaining fixed points and fiow properties that remained unexplored 
in [33j . In particular, we assess in detail the extent to which the renormalization flow of ^ inherits the structure 
of the perturbative KPZ flow. (Hi) Working on the reference to = 2 case for /i = 1, we prove explicitly that the 
additional interactions —Afk'^hk do not modify the asymptotic behavior for n = 3, 4; this justifies numerical results in 
[33} that exploited this fact and had remained analytically unproven thus far. As a bonus, this result also justifies our 
parameter restriction to the case /i < 2 in the definition of equation @ , since taking /i > 2 in the equation generates 
fc^/ik terms under renormalization, that become more relevant than the bare k^hi^ ones. 
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The paper is organized as follows: As a preliminary, we start in section |TT] with a scaling analysis of the family of 
equations (HI through power counting and through a more refined Flory-type approach \4^, and with considerations 
on their naive upper critical dimensions. Given the insufficiency of mere scaling analysis, in section ITTll we outline the 
perturbative DRG procedure as applied to equation We provide the proof on the mentioned vertex cancellation 
leading to the Galilean scaling relation. In section IIVI we discuss the most convenient parameter space in which to 
study the DRG flow. This is addressed in detail in section|Vl where we provide the main results (fixed points, stability 
and critical exponents) for dimensions d = 1, d = 2, and d > 2, assessing how dimensionality changes the overall 
properties of the equation. Finally, in section I VII we discuss the significance of our results both for the equations 
studied and from a more general perspective for surface kinetic roughening, and give our conclusions and an outlook. 
In order to make our presentation self-contained while trying not to obscure it, a number of technical details are 
provided in appendices. 

II. SCALING ANALYSIS AND UPPER CRITICAL DIMENSION 

In order to gain some intuition on the scaling behavior that can be expected for equation we can perform a 
scaling analysis. Straightforward power counting first leads to the conclusion that, under a rescaling of the height 
field and space and time coordinates as h{r,t) = b~°'h{br,b^) with 6 > 1, an equation with the same shape as (|4]) 
holds for h, with rescaled parameters 

D^b''-^'l^, JC^b'-"'JC, A = 6"+^^2^, flo ^ b'-^'^-'^Uo. (7) 

Here, a and z are, respectively, the roughness and dynamic exponents that characterize the scaling behavior of the 
system. The behavior of the ratio IC/i) immediately implies that /C is less relevant than ly under rescaling, since ^ < m. 
Likewise, under morphologically stable conditions, we expect a < 1, leading to the dominance of the non-local term 
over the KPZ term for /i < 1. Conversely, for morphologically unstable conditions we can assume a > 1, so that we 
would expect u to be irrelevant as compared with A for /i > 1. 

Power counting is indeed limited in its predictions regarding precise values of scaling exponents and expected 
behavior outside the exponent regions considered in the previous paragraph. One can improve over these results by 
performing a scaling argument a la Flory [49| in which one estimates the magnitude of individual terms present in 
(|4]). We assume that at long times t ^ t/, and averaged over length scales /, the height-height correlation function 
scales as C{l,t) ^ ft,^ ~ P", and that at long times the relaxation time of these fluctuations is ti ~ P. Note that ti 
and hi depend on the length scale I. The various terms in equation (jj]) can be estimated as 

z.(|j-ir/^k]|),^^jj, \{\{Vhf\)^^\^. (8) 

For white noise we can estimate its mean-square fluctuations on length scales I and times scales ti as rfi ^ {Hq/ SitiY^'^ , 
where Si is the average surface area of the interface on scale I, and thus Si ^ {hf + p)''/^ '49!], Yor morphologically 
stable conditions, the P term in 5*; dominates the height fluctuations (namely, a < 1); instead, for morphologically 
unstable surfaces the dominant term is h^. Let us start by considering the former case, and assume that v dominates 
over A. Then, by equating the non-local term with the inertial term (the one proportional to the time derivative of 
h), we obtain a simple condition for the characteristic time of the height fluctuations U ^ jv and, consequently, 
the relation z = [i. Under these conditions we estimate r\i ~ (Ho/Z'^ti)^/^, and equating in turn the inertial term 
with the noise fluctuations we obtain hi - l^i^-d)/^. This expression gives us the value of the roughness exponent 
a — (/i — (i)/2, which combined with z = /i leads to the so-called hyperscaling relation 2a + d = z, see section ITVl This 
set of exponents has indeed been observed in our numerical simulations of the morphologically stable condition [sS] . 
In the morphologically unstable case, we expect the noise term to scale as rji ~ {Ilo/hftiy/^ . If we try to equate this 
with the inertial term, we expect a violation of hyperscaling. Led by the previous power counting analysis, we now 
equate the inertial and KPZ terms, leading to hi ^ l^tf^ ^ ^ I", i.e. we obtain the so-called Galilean exponent 
relation a z = 2. 

In principle we have not assumed any restrictions on the values of /z when deriving the previous scaling relations 
and exponent values. As mentioned in the previous section, numerical simulations of the morphologically unstable 
condition 33] indicate that the scaling behavior depends non-trivially on the precise value of fi. Moreover, the present 
argument does not allow us to assess possible changes for different substrate dimension d. For these reasons, we need 
to resort to more refined scaling analysis, among which the DRG is a natural choice. 

Before leaving this section, we can also make estimates on the naive critical dimension of equation by making 
use of Eqs. ([7]). As in the standard KPZ case [s^l, the idea is to employ the exponents associated with the relevant 
linear theory in order to derive conditions under which A scales to zero, or else grows unboundedly for increasing b. 
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Typically such conditions depend on the value of d, hence the appearance of a critical dimension. In our case, given 
that the linearized equation ^ blows up exponentially in time for the morphologically unstable condition, such an 
estimate can only be meaningful, if at all, in the morphologically stable condition. Thus, employing the exponent 
values z = fi and a = (z — d)/2, we get A = 6"^A with nx — (3/i — 4 — d)/2. Hence, the nonlinearity becomes irrelevant 
for d > dciiJ.) = — 4, which is the case for d > 2 and any value of /i (note dc = for /i = 4/3). Alternatively, 
this implies that in d < 2 dimensions there exists a critical value /ic(c?) — (4 + (i)/3 such that for < /i < fic{d) 
the critical exponents are again given by the linear equation (mean-field). For fidd) < ^ < 2 the scaling behavior 
should be controlled by the nonlinearity; we reconsider this result in Section |V] below. In general, for > 2 we expect 
also KPZ-type behavior, but that remains beyond power counting analysis, as it requires renormalization. Note the 
strong similarities between these results and those expected for equilibrium systems with so-called weak long range 
interactions (36| , in particular about the possibility that non- localities can modify the values of the critical exponents 
(Slj . and the occurrence of mean- field type behavior as a function of /i. 



III. SKETCH OF THE RENORMALIZATION PROCEDURE 



Previous studies on the morphologically stable version of (|3]) (namely, > 0) for fJ. = 1 show that the equation 
indeed produces scale invariant interfaces whose critical exponents follow from simple dimensional analysis like the 
one performed in the previous section [7]. Actually, under conditions of morphological stability, one can even obtain 
analytically the exponents for arbitrary fi values in the deterministic limit that are expected to apply also for the 
stochastic case [30l|. However, the behavior for the morphologically unstable cases has been studied only recently [s^. 
In that work we integrated equation ^ numerically using a pseudospectral scheme 52, 53] for c? = 1, 2, and different 
values of fj,. The exponents measured numerically were confirmed through a DRG analysis using the Forster-Nelson- 
Stephen (FNS) [s^] scheme. In this section we sketch the general ideas of the latter approach, and show that the KPZ 
nonlinearity does not renormalize for any linear dispersion relation that is a linear combination of arbitrary powers 
of the wavenumber fc, a result that actually is independent on the morphological stability /instability of the system 
i.e., on the sign of v, and generalizes known results for the KPZ and Lai-Das Sarma- Villain equations. 

For the sake of generality, it is convenient to consider the evolution equation for a generic linear dispersion relation 
CTk. Thus, 

dM) - ak/ik(i) + ^■^[(V/i)2]k + %(i), (9) 

where the noise term is as in equation The DRG procedure starts by writing equation Q after time Fourier 
transform, 

[-(Tk - iuj] hk,^ = 77k,i^ -7:1 7^-\d [ q • (k - q) hq,n /ik-q.t^-n, (10) 

2 J|q|<A [^T^r J-oo 2n 

where w is the time frequency, A — vr/Ax is the wave- number cut-off in the system. Ax being the lattice spacing in 
real space. In Fourier space, the noise term in ([TU]) still has zero mean (77k, w) = and is delta correlated, but its 
variance is rescaled by a constant as 

{riK^ = 2Ho(2^)'^+i ,5k+k' (11) 

Following the standard FNS procedure [25l . [26l . Is^ l , the height and the noise fields are split into two types of com- 
ponents, fast modes h^^, Vkoj ^'^^ wave-numbers < fc < A/6, and slow modes h^^, Vuui ^'-"^ ^ k < A, where 
6 = e*' is a rescaling parameter (assumed to be larger that one, i.e. 61 > 0). Then one eliminates a "small number" 
(namely, 61 is infinitesimal) of the fast modes by solving the growth equation perturbatively for ^, substituting the 
solution into the equation for the slow modes and assuming statistical independence between high and low-frequency 
components. This procedure leads to an effective equation in which the fast modes are thus integrated out, 

[-cTk - S(k, 0) - i^] h<^ = _ ^ 1^ ^ I ^ q . (k - q) h<^ h<_^^^_,, + OiX% (12) 

where the effect of this coarse-graining step (the elimination of fast modes) is obtained by solving perturbatively the 
following integral [for the result in a particular representative case of equation see appendix \X\ 

I](k,c.)-2A2Ho^^ J ^[q-(k-q)][-q.k]Go(q,r!)Go(-q,-f^)Go(k-q,a;-r!). (13) 
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In Eqs. (IT2]) . ([T3)) we have denoted integrals over the fast (slow) modes with (/<)j and we have omitted the 

integration limits in the frequency domain [note that $7 S (—00,00)]. In equation (|13p. Go(k, cj) — [— (Tk — iw] ^ is 
the bare propagator, and the integral in q is calculated only for the (fast) modes within the shell q G [A/6, A]. In 
equation (I12p the coarse-grained propagator appears, 

G<(k,^) = [-ak-S(k,0)-ic^]-\ (14) 

From this expression, the effect of coarse-graining on the propagator is seen to merely amount to a modification of 
the system parameters. For example, E(k, 0) — Y.^k^ + Sk:^^ for equation ^ ior m — 2, so that the rcnormalizcd 
propagator has the same functional shape as the bare propagator, but with modified coefficients = + and 
/C^ = IC + Ek;- Actually, in equation (jl2p we have omitted higher order terms that can in principle be associated with 
a change in the value of A — >■ A< as induced by the present coarse-graining step, see section ITlI Al 
On the other hand, the coarse-graining of the noise variance is performed by starting from equation 

(/ik,>-k,-.) - 2n< G(k, c.)G(-k, -..), (15) 
in which the exact propagator appears. This equation leads, after a little algebra, to 

^k',.') = 2 [Ho + <i>(k, 0)] (2^)'^+! 5k+k' <5.+.' + 0(A3), (16) 
where the coarse-grained noise variance is given by 

^k,u;) = X'Ul J ^[q.(k-q)f|Go(q,17)|^|Go(k-q,c.-f))|2. (17) 

Hence, the noise variance in the equation for the slow modes is IIq — IIq -I- $. For the result in a particular 
representative case of equation Q , see appendix [Xj 

A. Vertex cancellation for generalized equations 

In this section we show that, to within our one- loop approximation, the coarse-graining procedure does not lead to a 
modification in the coefficient of the nonlinearity for equation Q with quite general CTk, namely, A^ = A. Specifically, 
we consider any linear dispersion relation that can be written as 

n 

a^ = -Y,c.\kr, (18) 

i=l 

where the set {r.i} is formed by a collection of positive real exponents that increase with the index, < r2 < • • • < r„, 
and coefficients ci are generic, but with the restriction that c„ is positive in order to obtain a well behaved dispersion 
relation at high wave numbers, i.e. CTk — > —00 for k — > 00. Given this choice, we consider an even more general 
equation of which ([9]) is a particular case, namely, 

at/ik(t) = <yMt) - h^T[iVh)\ + r^y,{t), (19) 

where ct > 0, see [11] for the special case in which ctj. — —vk^'^'^ . Note that equation ([T^ reduces to © for ct = 0. In 
the same way as the latter includes the KPZ equation as a particular case, equation (IT5)) generalizes also the celebrated 
Lai-Das Sarma- Villain equation [13, 13, that is obtained for /j, = a = 2. In short, the latter is a conserved version of 
the KPZ equation with non-conserved noise, and has played an important role in the study of kinetic roughening for 
thin films grown by MBE [2]. 

As is well known and is reviewed in detail elsewhere [5^ [2^ , the iterative solution of equations like (fTO|) or (fT9|) 
along the lines sketched above can be conveniently cast into a diagrammatic formulation (for notation and diagrams, 
see e.g. [13] )• Within such scheme, it is seen that the coarse-graining process is complete at one loop order once the 
coarse-grained coupling A^ of the non-linear term is evaluated from 




7 



where the original variables k and q in ^TU\i are replaced by the symmetric variables ki = k and k2 = q — k/2, and Tj 
are appropriate integrals that are evaluated next. The first one arises from the contraction of fast noise components 
with wave vectors ± + k2 — q) , and reads 

ri(ki, k2, .,,0.2) = iwfl^^' - ^ • - ^) + • - Y - 



'ki 



ki 
2 



k2 - qj Go (<?) Go (fci - g) Go U 



'j-kAcol'^ + h-q], (21) 



where we have used a shorthand notation, e.g., G ^fci — = G (ki — q, Oi — Q,). Using the change of variables 
j = q - ki/2, we get 



(g|ki-q|rq.(ki-q)(^+k, 



ki 



ki 



q-^-k2 — -k2 - — +k2-q 



ki 



J 



2a+4 



■ cos' (6ii) kf - cos' (62) + 0{kfkj), (22) 



where k^ • j = kj cos 9i. The contribution due to the bare propagator is a little involved, but we need to retain only 
the zero-th order contribution in ki of the perturbative expansion, which strongly simplifies the calculation. Thus, 



ki 



lim Go J + ^ Go - J Go J - fc2 Go fc2 - j 



+0(fci)+02 



where 



n 



^c,/0 +oik2) + n^ 



n-1 



Cnj' 



(23) 



(24) 



The function A^- evaluated at j = A takes a positive value A a due to the signs of c„ and a a {cn > and cta < 0, see 
above and section HV]) . Hence, after integration in 



+00 



(a| + z^)" 



2Af 



this contribution reduces to 

ri(ki,k2,o,o) 



dj 



•2(T+4 



(^_fc2) J> (27r)'^4A3(c„j'-")3 V4 



-cos' {61) k( - fc^cos' (62) 



2(f J> C^^r A3(c„j'-")3 7o 
A^Ho 



de sin"^' (6) [ - cos' (0) kf - ki cos' (6*) 



2d 



■(i+2(j+3 



(25) 



(26) 



where = 2/[(47r)'*/2p(^/2)] is the surface area of the c?-dimensional unit sphere Sd divided by (27r)'*, as defined in 
(|A8p . We stop the calculation at this point because the remaining contributions to the vertex renormalization in ([20)1 
will be seen to cancel (PSI) exactly. Indeed, the function Y2 is given by the sum of the terms obtained from the noise 
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contraction of ki — q with q — ki, i.e., 

2A2no f dq fdn 

r2(ki,k2,a;i,a;2) 



ki 



q ■ (ki - q) ( Y + ka 



q-y-k2 



(^^^-ka) -(q-ki) Go(g)Go (^i - g) Go ((Z-fci)Go L-'^-k 



(27) 



where 
9 



q-y-k2 



q-(ki-q)(^+k2 



ki 



ki 



q-^-k2 -^-k2 -(q-ki) 



•2(T+4 



and 



■ cos2 (6*1) kf - fc2 cos2 (6*2) + 0{kjkj). (28) 



(29) 



+ f j Go -ij Go (^i - Go (3- - ^2) = (c„,-)- [(A,^ + i^)(A, - ii)^ 
Following similar steps to the case of Fi, after integration in z, 



Az 



(A2 + z2)(A,±iz)^ 



4A3' 



we get 



A^Ho 



F2(ki,k2,0,0) = -^Xd / dj 



■d+2a+3 

A -3 



(30) 



(31) 



Finally, the contribution F3 in ()20p arises from a contraction of noise components with wave-vectors q and — q. This 
contribution reads 



r3(ki,k2,Wl,W2) 



2A2no 



f - kl\ J> (2- 



dq f dn 

2^ ll^^-^l 



^|.k2-q 



q • (ki - q) (-q) • ( y + k2 



ki 



ki 



qj Go (-?) Go (-g) Go (fci - g) Go 



fci 



+ k2- q 



(32) 



where 

|ki - q| 



ki 



k2 - q 



q • (ki - q) (-q) 



ki 



ki 



-k. 



ki 



k2 q 



( - cos^ (0i) kf - ki cos' (62) ) + Oikfk,), (33) 



and 



hm Go \ j 



+ y ) Go ( -3 - ^ ) Go ( j - I) Go (j - fc2) = (c„/")-' [(A| + z^)(A, + iz)^ 



ki 



(34) 



Thus, in parallel with the previous cases, in the wi,2 — >■ limit we obtain F3 = r2 = — Fi/2, hence from equation 
pO| finally A^ = A, so that the non-linear term does not renormalize at one-loop order. Notice, this result does not 
make any requirement on the morphological stability of the system described by equation , and has the analogous 
vertex cancellations for both the KPZ and the Lai-Das Sarma- Villain (LDV) equations as immediate corollaries. As 
is well known, for the former this property leads to the scaling relation a + z = 2 at the nonlinear fixed point, and is 
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usually associated with the explicit invariance of the equation under a Galilean transformation [U, [^iH^- However, 
for the LDV equation an analogous symmetry is known not to exist (soj . This fact makes it appear as a coincidence 
that the KPZ equation does combine such type of (Galilean) symmetry with vertex cancellation. Admittedly, vertex 
cancellation also occurs at two loop order in the KPZ equation [H^, [1^, but not in the LDV case fb^, so that the 
interplay between Galilean symmetry and the so-called Galilean scaling relation is subtle. Still, more recent works 
also question a causal relation between the two. Thus, it has been argued that for the Navier-Stokes (NS) and related 
equations, such as Burgers', Galilean invariance does not enforce vertex non-renormalization [soj. These arguments 
have been further substantiated more recently [60i . ,61.] within a field theoretical approach to the stochastic NS system. 
Likewise, recently models or discretization schemes of the KPZ equation have been shown to explicitly break Galilean 
invariance while fulfilling the a + z = 2 scaling relation, see [62l . Icsj and references therein. 



IV. DYNAMIC RENORMALIZATION GROUP FLOW AND FIXED POINTS 



Coming back to equation (U), in order to account for its scaling properties, we restrict our DRG analysis to the 
case TO = 2 and Af — 0, that is. 



(35) 



For clarity of the exposition, the details of the corresponding DRG calculations are left to appendix [X] Actually, 
the scaling properties of the equation do not change when we take into account a dispersion relation with additional 
stabilizing relaxation terms, less relevant than fc^, namely, n > 2 in equation (|4]). A proof of this is given in appendix 
|B]for the specific /i = 1 case, considering A/" ^ 0, with n = 3, and 4. 

After coarse graining of propagator and noise variance (as in appendix |A| , and the vertex (as in the previous 
section), the final step in the DRG method is the application of a rescaling that restores the initial value of the 
wave-vector cut-off. Moreover, this allows us to write parameter renormalization in a differential form, taking I as the 
independent variable. We can actually employ the formulae where the role of unsealed parameters is now taken 
by our coarse-grained estimates j/^, etc. and scaled quantities i>, etc. become finally the renormalized parameters. 
This makes explicit the non-trivial effect of the coarsening step as, e.g., differs from the bare parameter v in the 
original equation. After taking the SI ^ limit, the ensuing infinitesimal parameter variation constitutes the one 
loop RG fiow for arbitrary substrate dimension d (see appendix |A|) . 



dA . , 

- = - = A[a + z-2], 



"dT 



IC 



X^UoKd {d - 2)JC + {d- n)u 



Ad 



/C (/C + vf 



(36) 



(37) 



dllo 



2a -d- 



4 (/C + lyf 



(38) 



Without loss of generality, in equations (p6 | - (p8)) we have fixed the lattice cut-off, A = 1, bearing in mind that in 
order to avoid the finite pole of the bare propagator, the coarse graining step has been performed in an infinitesimal 
shell k E [A(l — SI), A] within the band of linearly stable (large) wave- vector values ^2^]. The dispersion relation for 
k = A (the largest wavenumber) can be written as 



CTA 



-ICA' 1 



K. 



-/CA^A. 



(39) 



where A^. = 1 + i^lC^^\k\^^^ is employed in the evaluation of the integrals in equations and p7p . The sign of cta 
gives us a criterion to identify a physically consistent DRG flow (see appendix for details). Note that cta has to be 
negative in order to have a well behaved dispersion relation and that, for the same reason, K. has to be positive. Thus, 
equation (j39p implies that Aa has to be a positive quantity. All the fixed points obtained within the DRG approach 
have to satisfy this requirement. 

Equations p6p - p8|) generalize the KPZ flow in a natural way for non-zero and inherit the known analytical 
limitations of the latter [li,[2|. Still, they carry valuable information. Introducing the following coupling parame- 
ters liil 



4(/C + ly) 



f = 



IC 



(40) 
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we can calculate their flow using equations ([36|) - (l38l) . Thus, 



dl 



(1 - /) {(m - 2)/ -^[(d- 2)/ + (rf - - /)] } , (41) 



^ = 3 {6/ - 4 - d + .g + 3^.(1 - /) + 3 I [(rf - 2)/ + (rf - ^)(1 - /)] } . (42) 

The ratio between the (possibly) morphologically destabilizing parameter, v, and the stabilizing one, /C, follows from 
the exact relation 1 — / = + JC). The condition Aa > is satisfied as long as these two couplings are positive 
or zero; for this reason we can restrict the DRG flow analysis to the first quadrant of the (/, g) plane. Then, a bare 
parameter choice corresponding to the morphologically unstable condition v{l = 0) < implies f{l = 0) > 1, while 
the converse condition f(l = 0) < 1 corresponds to morphological stability in terms of the bare parameters. The line 
f — 1 (at which = 0) is actually never crossed by any trajectory, since flow lines cannot cross each other in this 
autonomous flow and / = 1 is clearly invariant under the iteration. Note, / = 1 implies = 0, so that no explicit 
renormalization of v from negative to positive values is allowed by the flow (j4ip - (j42p . Thus, dynamical stabilization 
of the early-time morphological instability for equation (j35p and thus for equation Q cannot take place through 
Yakhot's mechanism (namely, "surface tension" renormalization) for the Kuramoto-Sivashinsky equation [24]. This 
agrees moreover with expectations based on a naive power counting analysis of both equations that would lead not 
to expect RG corrections to a fc'* term from a term that is irrelevant in comparison for /i < 2. 

The fixed points of the RG fiow are given by the values of the pair (/, 5) for which equations d//d/ — and 
dg/dl — hold simultaneously. Obviously, f — 1 and g = satisfy the first and the second equations, respectively. 
Moreover, the rhs of equation (|^T|) is identically zero for 

^ i^^-2)id-g)' ^^'^ 

while relation 

, ^ 3M./-l)+4 + rf-6/ 

leads to dg/dl = 0. 

A detailed analysis of equations (|4T|) - (|42|) yields the same fixed point structure as for the KPZ equation — namely, 
the Edwards- Wilkinson (EW) fixed points at (/,<?) = (1,0) and, for appropriate values of d (see next section), the 
KPZ fixed point at (/, g) = (1, 1/2) — , with the addition of two new nontrivial fixed points. One corresponds to the 
morphologically stable interfaces (we call it Smooth fixed point), being located at the origin of the {f,g) plane. Its 
dynamic exponent is z — fi, while hyper-scaling holds, namely, z — 2a + d. Note that a zero value for g means that 
the non-linear term A is vanishing and psp becomes a linear equation, the Smooth fixed point being the one employed 
in our discussion on the critical dimension in Section |lTl Actually, at the linear EW and Smooth fixed points, the 
equation becomes variational and hyperscaling ensues, relatiiig the dynamic and roughness exponents. Hyperscaling 
is associated with non-renormalization of the noise amplitude '26'] given that, under these conditions, the equation has 
(in the morphologically stable condition) the asymptotic height distributions V^wifA * cxp[(— /C/(2no)) / /c^|ft-kpdk] 
or T'smoothl/j] oc exp[(— i//(2no)) / fc''|ft,kpdk] at the EW and Smooth fixed points, respectively. 

The second new fixed point that exists for the fiow (|4T]) - (|42]) implements, as the KPZ fixed point, the scaling relation 
a + z = 2 usually associated with Galilean invariance, and we thus we call it Galilean fixed point. It is the one found 
in simulations for the morphologically unstable condition with v < [33l |. As in the Smooth fixed point, the dynamic 
exponent takes the value expected from power counting, z = fj,. The coordinates of the Galilean fixed point, (/*,5*), 
are given by equating (|43)) to (|44|) , namely, 

, -'^^t/^ y 9.-d + A-3^. (45) 
(^-2)(d-5*) 

Besides, this fixed point has to satisfy the conditions /* > and g^ > 0, hence it is not defined for some values of /i 
in the (0, 2] interval. Indeed, the condition > leads to an admissible value for this coupling for fi < (c? -I- 4)/3. 
For a fixed value of /i, this is equivalent to the d > dc condition using the value of the critical dimension derived in 
Section |TT] for the morphologically stable case. The second condition, i.e. /* > 0, has to be analyzed carefully. In fact, 
for d < 4/3, the Galilean fixed point occurs for /i within the intervals (0,d] IJ (4/3, {d + 4)/3] (for d—1, see the left 
panel of figure[T]). In the case 4/3 < d < 2, the admissibility intervals are (0,4/3) |J [d, {d + 4)/3] . Finally, for d > 2, 
we have only one condition for the existence of the Galilean fixed point, namely, fi < 4/3, as illustrated on the right 
panel of figure [TJ 
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In order to study the stability of these fixed points, we have to compute the derivatives of (|4T|) - (l42l) with respect 
to the couphngs, and evaluate the eigenvalues of the linear stability matrices. So far, we have outlined the main steps 
we have taken to obtain the fixed points of the DRG flow. In the next section we provide details on their stability for 
various substrate dimensions (in particular, d — I, d — 2, and d > 2). 

V. FIXED POINT STABILITY AND DEPENDENCE WITH DIMENSION 

The study of the trajectories of the RG flow provides useful information on the stability of the different fixed 
points. Although the stable case has been already examined in the context of the fractal KPZ (FKPZ) equation 
[3l| [w hich is obtained from ([4]) when v > and JC = Af = 0], our DRG study supplies additional information. 
In |3l] | the FKPZ equation was studied by means of a self-consistent expansion (SCE) for spatially correlated noise 
in the case of arbitrary substrate dimension. The limit of Gaussian white noise is recovered when the exponent of 
the noise correlations introduced in [sij is equal to zero. Although the SCE allows to identify the strong coupling 
KPZ fixed point with a reasonable accuracy (a = 0.295 and z — 1.705 for d = 2, compatible with the most reliable 
estimates for this value of d [fi^l), the SCE method has some drawbacks. Even though all the fixed points reported in 
the previous section appear within the SCE analysis, this method cannot assess which fixed point is selected by the 
dynamics during the growth process. For example, in the case of d = 2, three different phases are found [3l|, each 
one with different values for the critical exponents. For some value of the exponent of the fractional Laplacian (p in 
our notation, that is equivalent to 2 — p in [3l|), all these phases are admissible, there being no analytical criterion for 
the selection of one of them. The present section is precisely devoted to the analysis of the DRG predictions on the 
fixed point that is selected asymptotically by the dynamics for each substrate dimension, and as a function of bare 
parameter values. 

A. One-dimensional interfaces 

For d = 1 we obtain three fixed points whose coordinates do not depend on the exponent /i, i.e., the Smooth, EW, 
and KPZ fixed points (see left panel of figure [1]). On the contrary, the Galilean fixed point changes location in the 
{f,g) by changing the value of fi. In table |I] we summarize all the fixed points that exist for d = 1, their critical 
exponents, and the eigenvalues of their stability matrix. The (real parts of the) two eigenvalues of the Galilean fixed 

point, ^^^2' have been obtained numerically; we show in figure [5] how they change as functions of fj,. 



Name 


/ 


g 


z 


a 




Scaling relation 


Ai 


\2 


EW 
KPZ 

S 

G 


1 
1 


(5-3h)(1-m) 



1/2 


5 — 3/i 


2 
3/2 

M 


1/2 
1/2 

2 

2 - 


1/4 

1/3 

,,-1 

2m 
2-M 


2a + d^ z 
a + z = 2 

2a + d = z 
a + z = 2 


2 - ^ 
3/2-^ 
^ - 2 


1 

-1 

3/i — 5 


(m-2)(3m-4) 


M 



TABLE I: Fixed points in 1 + 1 growth dimensions (d = 1). In the first column, G stands for Galilean and S for Smooth, 
according to the terminology introduced in the main text. In the last column, Ai,2 are the eigenvalues of the linear stability 
matrix at the corresponding fixed point. Their real parts v'"^\, are plotted in figure [2] in the case of the Galilean fixed point. 
The latter is not defined in the intervals 1 < < 4/3 and 5/3 < fi < 2. 

In figure [3] we show eight representative snapshots of the RG flow for d = 1. For < /i < 1, the Galilean fixed 
point is unstable and moves on the sector of the (/, g) plane associated with morphologically stable bare conditions 
/ < 1, following the red line in the left panel of figure [T] Qualitatively, within this interval of /i, the flow behaves as in 
figures |3Ja)-(c), in which the DRG flow is shown for the cases fi = 1/4, 1/2, 1. Notice that the only finite fixed point 
for / < 1 and < n < 4/3 is the Smooth fixed point, as expected from power counting. For fi e (1, 4/3], the Galilean 
fixed point is not an admissible solution of the DRG fiow equations (pij) and (UU, and it disappears [for example, see 
figurelUd)] until fi becomes larger than 4/3, when it reenters the phase portrait as a fixed point for morphologically 
unstable bare conditions / > 1, see e.g. figure [3{e). Meanwhile, the basin of attraction of the Smooth fixed point 
changes as a function of ^, see figures [3lja)-(c). 
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For larger /i values in the interval 4/3 < < 3/2, the Galilean fixed point becomes stable, attracting all the 
trajectories with / > 1; this condition is shown in figure[31[e). The remaining trajectories are attracted by the Smooth 
fixed point. Actually, the Galilean fixed point merges with the KPZ fixed point at exactly /i = zkpz(c? ^ 1) = 3/2 
[see left panel of figured] and figure [HJf)], loosing stability in favor of the latter for larger values, 2:kpz(1) < A* < 2, 
justifying the observation of KPZ scaling in this range of /i, as obtained in numerical simulations under morphologically 
unstable parameter choices [33|. For 3/2 < < 5/3, the Galilean fixed point becomes a saddle, being located on the 
separatrix between the basins of attraction of the two stable fixed points for the DRG flow, the Smooth and the KPZ 
fixed points, see figure [3Ig), back on the morphologically stable sector. For /i > 5/3, the Galilean fixed point is no 
longer an admissible solution of the RG equations, and disappears. The only stable fixed point for these values of /i 
is the KPZ fixed point [see figure [3fh)]. As mentioned above, indeed it is KPZ scaling that is obtained in numerical 
simulations [33j performed for this range of /i values for morphologically unstable conditions. Actually, preliminary 
simulations performed for stable parameter conditions seem to agree with this conclusion also [65] , which in turn was 
expected from the power counting analysis performed in sectionllll above. In the first row of figure[4]we summarize the 
DRG stability of the four fixed point as function of /i. Note, numerical simulations for the morphologically unstable 
condition predict scaling as for the Galilean fixed point when < /i < 4/3 [ssl, and thus do not agree with the DRG 
fiow, probably due to the assumptions made in the perturbative expansion. Nevertheless, the DRG analysis gives 
significant hints about the critical behavior of this non-local class of growth equations, in particular the change of 
universality class from that associated with the Galilean fixed point for < /i < ZKPz{d — 1) to KPZ scaling for 
ZKPz{d — 1) < fj, < 2, generalizing in particular known results for the noisy Kuramoto-Sivashinsky equation [2Q. i21i]. 
that is the = 2 case of equation ^ . 

B. Two-dimensional interfaces 

The situation changes drastically for d — 2. As expected for perturbative approaches such as DRG [l^l and even 
field-theoretical schemes [stI. [ssj. there is no finite KPZ fixed point for the RG fiow, while the Galilean fixed point 
remains finite only up to = 4/3 < zkpz(2) — 1.61. For fj, = 4/3, this fixed point is at / = oo and g — 2, see the right 
panel of figured] As in the d = 1 case, the position of the Galilean fixed point depends on the value of fi, as shown 
in the right panel of figure dl In table |TT] we summarize the results for the fixed points in the two-dimensional case, 
showing in figure [H] the real parts of the eigenvalues of the Galilean fixed point, v{2j as functions of fi. We observe 
that the stability of the three fixed points does not change with fi (see the second row of figure [4|) . 



Name 


/ 


9 


z 


a 


13 


Scaling relation 


Xi 


\2 


EW 
S 
G 


1 


3{m-2) 





3(2 -m) 


2 




M-2 
2 

2-/1 




2m 

2-M 


2a + d^ z 
2a + d = z 
a + z = 2 


2 - ^ 
^ - 2 

(2) 




3 (m - 2) 

(2) 



TABLE II: Finite fixed points in 2 + 1 growth dimensions (d = 2). The meaning of G, S, Ai,2, and u'^' is as in table |T] The 
real parts wj^ji of the eigenvalues of the Galilean fixed point are plotted in figure (5] This fixed point is not defined for /i > 4/3. 

We can identify two main features in this DRG flow. The first one is represented in Fig. |6l^a)-(c). In this case the 
Galilean fixed point influences all trajectories in the morphologically unstable region / > 1; in fact, it is located on the 
manifold that separates two classes of trajectories: asymptotically, one type of trajectories flows towards f — I and 
g — >■ oo; the other type of trajectories flows towards / — > oo and g — 2, see flgures[6l^a)-(c). Within the morphologically 
stable region / < 1, the Smooth flxed point attracts some trajectories, while those outside the basin of attraction of 
this fixed point flow asymptotically to / = and g — > oo. 

For /i > 4/3, the Galilean flxed point is not an admissible solution of the flow equations. In this case all trajectories 
with / > 1 flow to / — > 1 and g — >■ oo, attracted by the KPZ flxed point at inflnity, see figure [Sl^d). The region 
with / < 1 does not display any qualitative change as compared to the case /i < 4/3. To compare with the present 
results of the DRG analysis, numerical simulations of the morphologically unstable case of equation ^ [s^ for d 2 
indicate Galilean scaling for /i < zkpz(2) — 1.61 and KPZ scaling for zkpz(2) < fi < 2. 
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C. The general case d > 2 



For d larger than two, we obtain a flow behavior that is independent on the substrate dimension. As for the KPZ 
case, there is no finite fixed point associated with KPZ scahng, while we do find an unstable fixed point located at 
/ = 1 (i.e., = 0) and g = gKr{d) that is associated with the non-equilibrium roughening transition (RT) [l[, hence 
our terminology. Starting from a bare parameter condition corresponding to /(O) — 1, the parameters of equation 
(|3]) flow towards the EW flxed point in case g{Q) < (7rT) or towards the strong coupling regime associated with KPZ 
scaling for g{0) > gjiT- See table HTll for the value of guT as a function of d. For this reason, the EW fixed point 
acquires a stable direction as compared with smaller values of d, changing from unstable (for d = 2) to saddle type. 
The two other fixed points. Smooth and Galilean, have the same stability and exponents as found for d — 2 (take into 
account that a negative roughening exponent means subleading corrections to scaling) . We summarize the results for 
d > 2 in table Iml 
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2a + d = z 
a + z ^2 
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TABLE III: Fixed points in d + 1 growth dimensions for d > 2. The meaning of G, S, Ai,2, and j is as in tables IT] [Til RT 
denotes the fixed point controlling the roughening transition as discussed in the main text. The coordinates of the Galilean 
fixed point depend on d as given by equations (|45l) . This fixed point is not defined for fj, = 3/2. 

As seen in the last two columns of table HiH the signs of Ai and A2 do not depend on d, implying the independence 
of the stability of these four fixed points on the substrate dimensionality (note that, for the RT fixed point, Ai > 
for every < /i < 2). The real parts of the eigenvalues for the Galilean fixed point are plotted in figure [7] for d — 3. 
As in the case of two-dimensional interfaces, one eigenvalue is positive whereas the other one is negative. We have 
checked that this behavior does not change by increasing the value of d. 

Beyond fixed point stability, and similarly to the KPZ equation, for any d > 2 the structure of the RG flow is 
almost identical to the two dimensional case, the only difference being that the / = 1 line acquires the RT fixed point, 
giving the possibility to observe EW scaling for z/ = 0. In figure |S] we report two representative flows for d = 3. By 
comparing them with those reported in figures [SJa) and [Sid), one cannot appreciate any qualitative change for flow 
trajectories outside the f — I line. Unfortunately, we are not aware of any numerical study of equation Q for d > 2, 
so that we cannot compare predictions of our DRG flow with numerical simulations for these values of dimensionality. 
In any case, we expect that the numerical results obtained for d = 1, 2 generalize for any substrate dimension. A 
qualitative summary of the DRG flows as obtained for different integer dimensions is provided in figure [SI 



VI. DISCUSSION AND CONCLUSIONS 



The DRG analysis of equation Q yields a number of non-trivial results that can be classified with respect to 
their relevance to several issues on surface kinetic roughening: (i) the fixed point structure for a class of non-local 
equations, and its dependence on dimensionality, as compared to available results from numerical simulations; (ii) 
general consequences for stochastic equations that feature a KPZ non-linearity, in particular the issue of Galilean 
invariance; (in) new mechanisms for the dynamical stabilization of morphologically unstable systems; (iv) the ensuing 
predictions on scaling behavior of realistic systems for which a continuum description can be expected to apply that 
fits within the framework of equation We will devote this section to a brief discussion of these issues, in a way 
that allows us to summarize the results of the paper and put them into a wider perspective. In general, interesting 
new results can be anticipated from further studies of the class on non-local equations ((H) along the lines discussed 
in the three subsections that follow. 
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A. Asymptotic properties 

For different dimensions, d, our DRG analysis can be roughly summarized as a generalization of both the KPZ 
and the noisy Kuramoto-Sivashinsky flows, in which additional non-trivial fixed points arise that are intrinsically 
associated with the non-local nature of the equations that we are studying. Thus, at the Smooth and the Galilean 
fixed points the dynamic exponent z takes a value that coincides with the one that would be expected from a naive 
power counting argument, while the roughness exponent can be evaluated from the latter through a scaling relation 
that is either hyperscaling or the one usually associated with Galilean invariance. 

The comparison between the analytical DRG results and those from numerical simulations leads to partial agree- 
ment. Specifically in the d = 1 case, the change in scaling behavior that takes place for increasing /i for the mor- 
phologically stable condition [6^1 is well described by the DRG results. Perhaps more remarkably, the DRG seems 
to also capture correctly the role of zkpz(c? = 1) —3/2 as a critical value for fi in the morphologically unstable case, 
that separates asymptotic behavior between that controlled by the Galilean fixed point for n < zkpz(1) and KPZ 
scaling for /x > 2:kpz(1), also in good agreement with previous numerical results [33|. However, possibly the most 
salient disagreement between the DRG and the numerics is the failure to generalize that role of Z}^pz{d) for general 
d > 1. This is not surprising in view of the known non-perturbative effects that are expected for the KPZ nonlinearity, 
particularly in higher dimensions. In general, the behavior of the Galilean fixed point changes with /x and system 
dimensionality in a rather non-trivial way, and some of the consequences for the DRG flow may be associated with ar- 
tifacts of the perturbative scheme employed. On the other hand, the dynamic renormalization group does account for 
the irrelevance of higher order contributions like those with parameter J\f in equation ([4]) , even in the morphologically 
unstable case, again as observed in numerical simulations [33l|. We expect that the same qualitative results apply to 
any equation with the shape of (jU using a general dispersion relation like ([T8|) such that ri < r2 = 2 < ■ ■ ■ < Vn- 

From a general perspective, the fact that equation Q — and more specifically equation — generalizes both the 
KPZ and the Lai-Das Sarma- Villain equations, has allowed us to prove vertex cancellation in the last two equations 
to stem merely from the shape of the nonlinearity, and to be independent of an underlying continuum symmetry of 
the dynamical equation. Although limited to the one-loop perturbative expansion that has been employed in our 
DRG approach, this result may prove informative in the current questioning of Galilean invariance as inducing vertex 
non-renormalization in the KPZ equation, although further evidence is needed. In general, a more detailed analytical 
understanding of the asymptotic properties of equation (j4]) will require resorting to techniques that differ in nature 
from the perturbative DRG. 

B. A mechanism for dynamical stabilization 

A very remarkable feature of the present DRG results is that they describe analytically a mechanism for the 
dynamical stabilization of morphologically unstable systems, that differs from Yakhot's stabilization mechanism for 
the (noisy) Kuramoto-Sivashinsky equation. Namely, the latter requires renormalization of the negative "surface 
tension" parameter z/ to a positive value, eventually leading to asymptotic behavior in the KPZ universality class. 
Dynamically, as seen in numerical simulations of the noisy Kuramoto-Sivashinsky equation both in one [l9, 21, 6^ 
and in two |22| dimensions, the periodic cell structure that appears at short times is seen to evolve in a chaotic fashion 
that leads to height disorder and kinetic roughening properties at sufficiently long times. In marked contrast, for 
equation Q the corresponding unstable v parameter cannot change sign for /it < 2, as implied by the first equation in 
(1551) . or equivalently by the fact that / = 1 is an invariant of the DRG flow, see equations (|in)) - (|^T|) . Still, the system 
evolves from an initial periodic morphology into a late time topography that is also disordered and rough. For any 
value of /i < 2, this process is more similar to cell coarsening by creation and annihilation of cusps than to the familiar 
chaotic cell dynamics characteristic of the (/i = 2) nKS case. See movies, e.g., in the Supplemental Material provided 
with (43l] and with [331] for one and two-dimensional systems, respectively. Cusp annihilation and creation has been 
studied in detail in the particular case of the stochastic and deterministic Michelson-Sivashinsky equations, namely, 
for /i = 1 with m = 2, see e.g. [G^I and references therein, but does not seem to have been studied explicitly for other 
values of fi. Note that, in spite of such strong differences in the dynamical evolution between the nKS equation and 
equation Q, they share the same KPZ asymptotic state, provided v < and (as conjectured) zkpz(c?) < /i < 2. 

C. Comparison with experimental and other model systems 

In view of the previous results on scaling properties of equation ^ , a natural question concerns their experimental 
relevance. We believe that the kinetic roughening properties of a number of experimental systems are indeed com- 
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patible with the predictions from equation For instance, in the morphologically stable condition, viscous flow 
induced on glassy films by ion irradiation leads to the scaling exponents obtained from @ ior fi — 1 |68|] . 

Perhaps more interestingly, a number of actual thin film systems that are characterized by large values of the growth 
exponent f3 = a/z [l| are compatible with our continuum framework in the morphologically unstable case, possibly 
due to the occurrence of geometrical shadowing effects 40]. For instance, in 69] Si(lOO) surfaces were plasma etched, 
finding a — 0.96 ± 0.06, (3 = 0.91 ± 0.03, and, z — 1.05 ± 0.09, consistent with the predictions of equation ^ for 
/i = 1 at the Galilean fixed point. Very similar exponent values have been recently found in chemical vapor deposition 
(CVD) growth of silica films p^ . CVD is an experimental technique within which shadowing instabilities are quite 
generically expected. Thus, e.g. in [tO] amorphous hydrogenated silicon (a-Si:H) was deposited on glass via a plasma 
enhanced CVD technique, finding /? = 1.13 ± 0.04 for the condition of low-silane concentration in the gas mixture. 
Similarly, (7l| reports the roughness evolution of hydrogenated amorphous carbon (a-C:H) films grown by electron 
cyclotron resonance CVD. A growth exponent value /3 = 0.74 was measured, which suggests that an underlying 
morphological instability is responsible for the cauliflower-like morph olog y observed at the late stages of the growth 
process, that looks quite similar to that found in Further, in j72l | a large growth exponent (/3 = 0.9 ± 0.1) 

is reported for hydrogenated diamond-like carbon (a-C:H) coatings deposited through plasma enhanced CVD. In 
general, it would be interesting to search for experimental parameters/techniques that allow to tune the value of fi 
characterizing a putative continuum description through equation ([4]) , in order to probe the corresponding continuum 
line of fixed points. A potential candidate might be the sticking probability in CVD systems, that seems able to tune 
scaling properties from KPZ-type [tI] ("large" /i values) to sMS-type (/x = 1) 0, although values of the effective ^ 
parameter should be still identified that differ from fi — 1,2. 

In connection with kinetic roughening in the presence of shadowing instabilities, we note that other continuum 
and discrete models are available [3] . Typically many of them involve the calculation of the so-called exposure angle, 
which is a non-linear and non-local function of the height profile, that is very expensive to compute. By means of a 
small slopes expansion, the exposure angle can be seen to lead to a term of the form khk in the height equation [t^ I , 
leading to a non-local dispersion relation as in (U)) (with < 0, /i = 1, and, to = 4). Note that for this choice of /i, the 
Galilean fixed point leads to the exponents a = l3 — z = 1. However, other authors have included the non- locality of 
the shadowing process in the non-linear part of the evolution equation through a non-local redistribution process of 
the incident particle flux [ill . The latter model has been studied numerically with a Monte Carlo method [tI, [t^ , 
and analytically (77j . finding the same critical exponents, i.e., a ~ /? ~ z ~ 1, and a = 1.04, /3 = 1.08, and, z = 0.96, 
respectively. In spite of the different formulation of the shadowing instability, the critical exponents coincide. 

Overall, we believe that equation (jH) provides a very compact description of a very wide class of systems undergoing 
non-conserved interface growth, both under conditions of morphological stability and in the presence of pattern forming 
effects. Note that, in spite of the presence of the KPZ nonlinearity in the equation of motion, KPZ scaling only occurs 
under the, relatively restricted, simultaneous conditions of morphological instability and (as conjectured) sufhciently 
large /j, > ZKPzid) values. This fact, together with the long transients induced by morphological instabilities, may 
account for the surprising experimental scarcity of KPZ scaling Understanding the degree to which this is truly 

the case will require further exploration of the present and related continuum models. 

Acknowledgments 

We gratefully acknowledge discussions with C. Escudero. This work has been partially supported through Grants No. 
FIS2009-12964-C05-01 and No. FIS2009-12964-C05-03 (MICINN, Spain). M. N. acknowledges support by Fundacion 
Carlos III (Spain) and by Fondazione Angelo Delia Riccia (Italy). 

Appendix A: Coarse-graining of propagator and noise variance for Eq. (|35[) 

For simplicity and as justified in more detail in the next Appendix, we consider equation (|35|) . namely, equation (|4]) 
with (Tk = —vk^ — K.k'^ , as the simplest representative for the class of non-local equations studied in the main text. 
Here we present the detailed calculations for the coarse-graining step in the renormalization of the propagator and 
the noise variance for equation p5p . 

In equation (|12p for the slow modes, the contribution from the elimination of the fast modes is contained in the 
function E. In order to evaluate expression (1131) . we symmetrize the wave- vector integral by introducing the standard 
substitution (q, fi) — >• (j + k/2,ri -I- uj/2) [53|. We begin the calculation by considering the various contributions 
separately. 
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The wave-vector contribution due to the vertex is readily obtained in the new symmetrized variables as 

[q • (k - q)][-q • k] = f cos(0)fc + ^fk^ + 0(fc3), 



(Al) 



where 9 is the angle between k and j. Next, we expand the contribution from the bare propagators by taking the long 
time limit (i.e., w — >■ 0), 



hm lGo m'Go (k-q) 



2\ m/2 



uj'' [1 + x cos{9) + + ICf (l + X cos{e) + ^ 



,2 \ m/2 



2\ 



vj^" ( 1 - xcos(6') + ^ j + /Cj2 M - a;cos(6l) + ^ 



2\ M/a 



Ai +iz Aj + z 



2A, \^ ,^,k 



{K:ff{A] + z') (A,+iz) 



(A2) 



where we have introduced variables z — Vt/Kp and x — k/j, and the function Cj — 1 + i^/i/2/C|jp^^ in order to 
shorten the notation [A^ has been already defined in ([M)) ]. We can multiply the last result by (|Aip and split the 
integral in q into an angular integral in d— 1-dimensions and a one dimensional radial integral, where we approximate 
the shell to a (hyper-)spherical shape. Thus, 



J> (27:)" (27:fJ> Jo 



l> (2^)'^ (2^)'' J> 

where Sd = 2it'^^'^ /r{d/2) is the surface area of the d-dimensional unit sphere. 
We can rewrite (|13p up to second order in the external wave- vector fc. 



(A3) 



E(k,0) = 2 



•^^IIo Sd-i 



K. (27r 



,d+l 



A/6 



[ de* sin'*~2(g|) 

-'0 



J cos(0)/nfc + ( ^/ii + C, cos2(0) (/i2 - 2A,/2i) ) ^ 



where the integrals 



lu 



+ 00 



dz(A2+^2N '(Aj.+i^)- 



(A4) 



(A5) 



are easily calculated through the residue theorem as In — 7r/2A|, /12 = 7r/4A^, and I21 = Btt/SA^ (taking into 
account that A^ > for j € [A/6, A]). 
Integrating (jA4p in 9 using 



/ d6'sin''-2(6i)cos(6') =0, 
"'0 

r <19 sin'^-^{9)cos''{9) = \ C &9 sin'^-^{9), 
Jo dJo 



we obtain 



Finally, we perform a perturbative expansion up to first order in 1/6 = e^"^' = 1 — SI + 0{SP), that leads to 



E(k,0) = 



4d (/CA2-M + 1,) 



[{d - 2)K.^^-^' + {d- 61 k^ 



(A6) 
(A7) 
(A8) 

(A9) 
(AlO) 
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As for the coarse-graining of the noise variance, we proceed as above using the symmetrized variables, 



lim |Go (g)r 



Go(k-q 



[Aj + Cj cos{9) X + z^] [Aj - cos{0) x + z^] \ 



Putting this last result into the expansion of equation (jl7p up to zero order in powers of fc, we obtain 



$(k,0) 



d-l 



+00 



dz 



Sd-1 



^ICff (Aj+z^y Jo 



4/C3 



A/b 



A3 



SI. 



4 (/CA2-M + ,yf 

The modified parameters are obtained from equation (jA10[) and equation (jA12[) 
v< ^v, A< = A, 

A^Hoi^d A'^+^-^A' 



Ad (/CA2-A' + vf 



[(d - 2)/CA2-'^ + {d- 81 



n< = Ho <^ 1 + A^Ho 



4(A:A2-M + iy)- 



(All) 



(A12) 



(A13) 



Now we can perform the rescaling according to the usual DRG transformation. Proceeding as indicated in the main 
text, the parameter flow finally reads 

di^ d\ 

- = - = A[a + z-2], 



d/C 
dT ' 

dllo 



lClz-2 



Iio{z-2a-d + A^Ho 



4d (/CA2-M + vf 



[{d-2)K.K^-^' + {d- ^)v] 



4 (/CA2-M + jy)- 



(AM) 



The set of equations (f36|)-([38l) in the main text follows from (jAl4|) by imposing an unit wave-vector cut-off, i.e. A = 1. 



Appendix B: Irrelevant stabilizing terms 

In this appendix we present the DRG calculations for a case of equation ^ in which the linear dispersion relation 
includes additional stabilizing terms (proportional to and fc"*) for 1-1-1 growth dimensions. We will see that the 
linear terms parametrized by A/" ^ and exponent n larger than 2 are irrelevant with respect to a term like —k^h\^, and 
do not change the hydrodynamic scaling behavior of the equation. This result can be generalized straightforwardly 
to the full family of nonlocal equations considered in the main text, although for the sake of concreteness we provide 
the detailed DRG calculations for the case /i = 1 only. This result is not surprising, since these terms are linear, and 
simple dimensional analysis anticipates such results. However, we have included these details for completeness, and 
to explicitly show how the non-trivial fixed points that control the dynamics of equation Q already emerge when 
considering only the "surface tension" term proportional to k"^ , as done in sections |IV] and El and appendix |^ 

To begin with, we consider a linear dispersion relation of the form 

<Jk = ~Ci|A:| - C2k^ - Cslfcp - C4fc^ (Bl) 

so that in the coarse-grained propagator we have to retain powers of wave-vector up to fourth order (note that here 
k is the one-dimensional wave- vector and not its modulus). Later on we will consider the choices C3 = and 6*4 = 
as particular cases. We proceed as in appendix O using the standard symmetrization and considering one by one the 
various contributions to the integrand in equation (1131) . The first contribution is readily obtained as 

k{q - k)q' = k^j + ^) - ^) = fk + - - \k\ (B2) 
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The product of the first two bare propagator factors in equation ((T3l) is evaluated considering that |fc| <§i \ 

b- + fc/2|-b1+sgn(j)A;/2, (B3) 

|j + fc/2|3 = |jf + 3j2sgn(j) fc/2 + 3\j\kyA + sgn{j) e/8. (B4) 

Using these identities, we are ready to expand the first product of two bare propagators appearing in up to third 
order 



lim IGq (9)1^ = — 7 

-^0 {Ciff 



1 



(Aj + ajX + fej^E^ + CjX^^'^ + 



1 



1 



A? +z2 



A2 + z2 



(2ajAj)2 



2fe, Aj - a2 



A2 + z2 



(26,A, + 4) _ 2(c, A, + a,b,) - / /" ^"^ 

(A| + z^j 



(B5) 



where 



A, = l 



C2 



3 C'a C2 Ci 



2 QUI ' C4j2 ■ 2C4bf' 



3 3 C. 



C-2 



2^ ACM ACiP' 
1 , C3 



2 8C4|jr 



(B6) 
(B7) 
(B8) 
(B9) 



and we have introduced two new variables: the non-dimensional expansion parameter x = k/ j, and the rescaled 
frequency variable z = 0,/ (C4J'*). The same procedure as applied to the remaining factor in the integrand of equation 
(ini) leads to 



lim Gn ( k 

^0 V 



77 (Aj — fljX + hjX^ ~ CjX^ + iz) 



G4j4 

1 1 

C4 Aj + iz ' 



Aj + iz 



A, + iz \ A, 



b, x^ 



A, + ii 



2ajfej 



Aj + iz (Aj- + iz)^ 



(BIO) 



After multiplication of (jBSP and (|B10p . we can integrate out the rescaled frequency after the change of variable 
O — > z = ri/(C4j^). Note that the frequency integral in the new variable is 



+00 



dl7 ^ C^f 



+ CXD 



dz. 



(Bll) 



We can use the compact notation introduced before for the integrals involved in the evaluation of the coarse-grained 
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propagator [see equation (jASp j. After some algebra, we obtain the contribution to equation (IT3|) as 



+00 



1 



+ 



Cj/12 - 2aj6j/i3 + a|/i4 - 2aj Aj (0^/23 - ^^^22) + (2ajAj)^ 132 - 2bjAjl22 - a^jl. 



+ 4aj Aj (26jAj + a^) /31 - 2 (c^ A^ + a^bj) hi - {ia-jA-jf hi 



(B12) 



The Iim integrals are easily calculated using the residue theorem, giving the contributions listed below: 
hi = 7r/2A2, /12 = tt/AAI /i3 = tt/SA^, hi = Tr/IGA^, hi - Stt/SA^, 
/22 = ^/4A5, /23 = 5^/32A6, /31 = 57r/16A6, 732 = 15^/64Aj, /41 = 357r/128A8. 

Using these relation and multiplying contributions (jB2l) . (|B12p together, we can write S as sum of four terms. 



where 



fiU) 
Uij) 



E(fc,0) 



(CiA.ff 



1 

{CiA.ff 
1 



47r 



2 A, 



1 

A" 



2A. 



/2(j) = 




1 ^ 


(C4A,j3)2 V 




26,^ 


a, 
2A, 


(aj 
I A. 


-26, 


]+^ 
1 4A, 



.A . Q I ^ a/^ ^ A2 



(B13) 

(B14) 
(B15) 
(B16) 



Considering the parity of the various functions of j that appear, we can simplify expression (|B13p . In fact, aj, bj, Cj, 
and Aj in (|B14l) - (jB16p are even in j: for this reason the functions /i, /s are odd, whereas /2 and fi are also even. 
The momentum shell is a symmetric interval around j equal to zero, hence the integrals of fi{j) and fs^j) vanish. 

Another important point to address is the contribution of the shell to the renormalization of the propagator. This 
calculation is only required when k'^ (or higher order) terms are present, because the lowest contribution from the 
shell occurs, precisely, at that order. Writing explicitly the integration limits in the positive part of the momentum 
integral, and expanding for small k [t^, [t^ , 



A-k/2 



A+k/2 



J djf{j) = ^ J dj7(j)+ J djfif) 



A/b+k/2 



A/b-k/2 



A/b 



djfU) 



-/ A 



16 



/' A 



A k 
5 + 2 



/ A 



2 ■' \b 2 



, + /' A + - 
2 ■' V 2 



A k 



(B17) 



fc=0 



where / is any function of j (the symbol + in the integral means that the integration interval is restricted to positive 
j). With this simple result, we see that if a function /; in equation (jB13p contributes at order n in powers of fc, then 
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through the expansion of the shell it also contributes at every order n + 2m, with m a positive integer. Thus, in our 
calculation we only need to consider the contribution of /2 to the fc^ term as induced by the shell. 



A/b » 



/2(A) + -/^'(A) 



ASl. 



(B18) 



Here, the function F2 is the primitive of /2. We now introduce three new coupling variables that help us to write 
(jBlSp in a more compact way, 



F = 



B = 



so that we can express /2, its second derivative, and using these variables, 
1 3 + 2B + F 



f2{j) 



2{CifY {l + A + B + F) 
1 



[11 + A{)A + 27 B + 29F + 2A^ + 195^ + 11^^ + A7AB + 21 AF + S5BF 
+ {2B + F) {A^ + 2^2 + + TAB + iAF + ABF) [1 + A + B + F)~' 



(B19) 



(B20) 



(B21) 



/2'(J) 



■ [63 - 90A + 119B + 27F + 9A^ + 82B^ + 15F^ - 95AB - 27 AF + 55BF + 3F^ + 2A^B 



3AF^ + 20B^ + UBF^ - 26AB^ - 17 ABF + 25B^F] {l + A + B + F) 



-5 



(B22) 



Considering that the shell is composed of positive and negative wave-vectors we have to count twice the previous 
contribution. 



E(fc,0) 



An 



2/2(A)fc2+(2/4(A) + i/^'(A))fc 



ASl = (E2fc^ + ^ik"^) SI, 



(B23) 



where and S4 are 



A^Bo 



(3 + 2B + F) 



S4 = - 



47r C|A5 (l + A + B + F) 
A^Bo 



3 ' 



167rC|A7 



52 - 130A + 92B - 2F + + 63B^ + AF^ + 16B^ + 2F^ - 1A2AB - ASAF 



(B24) 



20BF - 3QABF ~ AOAB^ - GAF"^ + 8BF^ - A^F + IbB^F 



(l + A + B + FY 



(B25) 



In our DRG program, the next step is calculating the contributions arising from the coarse-graining of the noise 
variance. We start by considering equation (IT7|) for d = 1. On the other hand. 



(B26) 



while the first factor due to the bare propagator has been already obtained, see equation (jB5|) . The other contribution 
is exactly the square of the absolute value of equation (jB17p . 



lim 



Go 



2ajAj 
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(B27) 
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The coarse-graining of the non-conserved noise involves only the zeroth order term of this expansion, hence 



(B28) 



After integration in z and, due to the parity of the functions in (|B28p . we can finally write the coarse-grained noise 
variance in the following form, 



$(A:,0) 



47rC| J A/6 



A^n^ SI 



j8A3 47rC| A7A3 ' 



(B29) 



As seen in the main text, there is no renormalization for the KPZ vertex for any dispersion relation that can be 
written as a polynomial in fc, at least within our one- loop expansion. Now we are in a position to apply the usual 
rescaling, and to calculate the renormalized parameters 
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1 + A^Ho 



1 + A^Ho 
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47rC2C|A5 
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where the functions P2 and P4 are 
3 + 2B + F 
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(B30) 
(B31) 

(B32) 

(B33) 

(B34) 
(B35) 



In the next two sections we show that for a dispersion relation as the one appearing in equation (U) with m = 2, 
additional relaxation terms with n > 2 do not change the hydrodynamical properties of the system as obtained in 
appendix El 



Irrelevance of the k term 



As a first case we consider the dispersion relation that is a third order polynomial in | fc | , (7^ = — Ci | fc | — C2 fc^ — C3 1 fc ^ . 
The ensuing evolution equation is an interpolation of the sMS and MSKPZ equations. Using the results already 
obtained, we can write the RG flow. 
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(B36) 
(B37) 

(B38) 
(B39) 



From equation (|B37p . note that, even if the bare "surface tension" coefhcient C2{1 = 0) = 0, a non-zero value for C2 
will be generated by the RG flow, provided Cz{l = 0) 7^ 0. Introducing the coupling variables 



ai 



C3A2' 



C2 

C3A' 



9i 
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the resulting flow reads 
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dl 
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/i(l + ai + /i) 

91 

{l + a, + hf 



In this case the condition Aa > that we have used in order to calculate the integrals requires 1 + ai 
In fact, when the dispersion relation is a third order polynomial, the coefhcient C3 has to be positive and 



CTA 



-C3A3(l + ai + /i), 



(B41) 
(B42) 

(B43) 
/i > 0. 
(B44) 



has to be negative in order to have a well-defined short distance behavior, hence the positive sign required for the 
parenthesis in (jB44l) . For this reason we do not consider fixed points of equations (|B41[) - (jB43p that do not satisfy this 
restriction. The only admissible fixed point in the (ai,/i, gi) parameter space, Fpa^Oj is reported in table ITVl This 
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ai 


fi 


91 
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Ai 
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1 


1/3 


2 


1 


4 



TABLE IV: Fixed point of ((B4T1) - (|B43) ) in the (ai, plane. 



fixed point is associated with the linear limit in which Ci = C2 = A = 0. All the eigenvalues of the stability matrix 
associated to Fpa^o are positive, hence this fixed point is linearly unstable in every direction. 
A second set of couplings that can be considered is 
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their flow reading 
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(B45) 

(B46) 
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(B48) 



The corresponding fixed points and the eigenvalues of their stability matrices are listed in table |Vl whereby they are 
seen to correspond to the EW (Ci C3 = A = 0) and KPZ (Ci = C3 = 0) fixed points. Note they both have stable 
and unstable directions. 



Name 


a2 


h 


52 


z 


Q 


P 


Ai 


A2 


A3 


EW 











2 


1/2 


1/4 


1 


-1 


1 


KPZ 








1/2 


3/2 


1/2 


1/3 


1/2 


-3/2 


-1 



TABLE V: Fixed point of ((B46)) - (|B48l ) in the (02, {2,92) plane. 
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A third convenient set of couplings is 



0.3 



their flow equations reading 



das 
dl 

dh 
dl 

dga 
dl 



-2a3, 

ff3(2a3 + /a) 
(1 + as + hf 

93 



53 



/3, 



53 



(1 + ^3 + /3)' 



(B49) 

(B50) 
(B51) 

(B52) 



In this parameter space, the Smooth and the Gahlean fixed points that have been discussed in the main text arise 
naturally, with features as in table IVll Note, the Smooth fixed point is completely stable in every direction. 



Name 


^3 


h 


93 


z 


a 


13 


Ai 


Aa 


As 


Smooth 











1 








-2 


-1 


-2 


Galilean 








2 


1 


1 


1 


-2 


1 


2 



TABLE VI: Fixed point of ((B50l) - (|B52| ) in the {as, /s.ffs) plane. 

Using the values for exponents a and z as calculated for each fixed point, we can show the irrelevance of C3 
compared to C2. After rescaling, equation Q reads 



as 



dthk 



-Ci\k\ 



2^2 



C3 
62 



A 



J- [(V/i)2] +6(1 



-2a)/2 



(B53) 



and we consider now the case of C2,C3 > 0. Clearly, from linear analysis the ratio between the stabilizing terms 
is C26/C3, so that the surface tension term becomes more relevant than the |fcp term. Moreover, the sign of the 
eigenvalues of the linear stability matrix at the EW fixed point shows that this fixed point is stable along the 
direction /2 — C3A/C2. The same situation occurs for the KPZ fixed point (Ci — 0), in fact it has two negative 
eigenvalues associated to the two stable directions /2 and 52- For the Galilean fixed point, the only stable direction 
is fli — CsA'^/Ci, while for the Smooth fixed point every direction is stable. The linear stability of these fixed points 
shows that Fp3^o can be reached only if Ci = C2 = A = and is unstable in every direction, hence we can conclude 
that the term is irrelevant as compared to (for C2, C3 > 0). We would need this term only in case of a negative 
surface tension, i.e., C2 < 0, in order to have a dispersion relation that is well posed in the — >■ (X) limit. 



Irrelevance of the k term 



In the last part of this appendix we show that the parameter C4 is irrelevant when the dispersion relation con- 
tains a relaxation term with an exponent smaller than four. Let us consider a dispersion relation of the form 
CTfc — — Ci|fc| — C2k^ — Cik'^: in this case, expressions (jB24[) and (|B25[) reduce to 



S2 = 



A^no 

' 167rC|A7 



F 



{l + A + F)-_ 
52 - 130A + ~2F + 2F^ 



A8AF + 7A^ - 6AF^ - FA' 



{i + A + Fy 



(B54) 
(B55) 
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Hence, after rescaling, 



dCi 
dl 

dC4 
"dT 



dllo 
dl 



= Ci [z - 1] . 



dA 
dT 



= X\a + z-2] 



C4 



z-2 



47rC2 (Ci +C2A + C4A3)^_ 



(B56) 
(B57) 



z-A 



(52C|A« - 130CiC|A5 - 2C2C|A*^ + 7C^C4A^ + ACldM^ 



167rC4 

2C|A2 - 48C1C2C4A3 - 6CiC|A - ClC2^ (Ci + C2A + C4A3) 
A^HoA^ 



z - 2a - 1 + 



47r 



(Ci + C2A + C4A3) ' 



Analogously to the result obtained in the previous section, equation (|B57p implies that, even if the bare 
tension" coefficient C2{1 = 0) = 0, a non-zero value will be generated for C2 by the RG flow, provided Ci{l = 
Also in parallel to our previous computation, we consider a first set of couplings 

_ Ci f _ ^2 _ A^IIo 

for which the renormalization flow is 

3 - ^ (52 - ISOai + 4/2 - 2/1 + 2/3 - 48ai/i + 7a? - 6aiff - /la?) (1 + ai + /i)~ 



dfli 
"dT 



ai 



d/i 
d^ 



2/1 + ^ (12 - 24/1 + 22/2 _ 2/4 + 24ai + 12a? + 56ai/2 + 162ai/i - 3/ia? + a?/^ + Gai/^) 



X (l + ai+/i)-% 



(B58) 

(B59) 

'surface 
0)^0. 

(B60) 
(B61) 

(B62) 



^1 
dl 



91 



(i+«i + /ir 



1 - ^ (52 - laOai + 4/2 - 2/1 + 2/f - 48ai/i + 7a2 - Gai/f - f^af 



X (1 + ai + /i)" 



(B63) 



The different flxed points in the (04, /i,5i) parameter space are given in table IVlIl where the coordinates of the KPZ 



Name 


fli 


/i 


91 


2 


a 


/3 


Fp4,o 











4 


3/2 


3/8 


KPZ 





13.1868 


1064.4 


1.54 


0.46 


0.299 


FP4,1 


7.17 


-6 


20.4 


1 


1 


1 



TABLE VII: Fixed points of (fB6T|) - (|B63l ) in the (ai, plane. 

fixed point have been determined using the only admissible real solution, xq, of equation x^ — 8a;'^ — 63a;2 — 68a; — 42 = 0. 
Hence, the KPZ fixed point is located at the point (0, xo, P(xo)), where P{x) = (686 + 969a; + 846a;2 - 25a;3)/97. 
Likewise, the Fp4_i fixed point has been obtained solving 7y'^ + 2?/ — 374 = and using the only admissible solution yo 
for the components of vector {yo, —6, Q{yo)), where Q{y) = (13014y — 92286)/49. The Fp4_o fixed point is associated 
with the linear Molecular Beam Epitaxy (MBE) equation (Ci = C2 = A = 0) [H, while Fp4^i has the same exponents 
as (but does not coincide with) the Galilean fixed point. The stability matrices of these points are dense, and we 
report their eigenvalues in table IVHII All these points are unstable, or at least hyperbolic, in the three dimensional 
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Name 


Ai 




As 


Stability 


Fp4,0 


3 


2 


7 


Unstable 


KPZ 


0.54 


-2.90 


-1.06 


Saddle 


Fp4,l 


88.4 


7.33 


-0.24 


Saddle 



TABLE VIII: Eigenvalues of the fixed points of table lylll 



space, but if we consider the Ci = plane, the KPZ fixed point is stable. 
A second set of couplings is provided by 



02 



Cj_ 

CoA' 



f2 = 



C4A2 
C2 



whose flow is 



da2 
"dT 

dl 



= 0.2 



= h 



1 - 



92 

4/2 



52(1 + 3/2) 



dg2 
dl 



92 



(1 + 02 + /2)3 

2 - 12/1 + 24/1 - 162a2/| - 22/| - 5602/2 - 602 - aj - Uajf^ - 24a2/| 

X (1 + 02 + /2)" 
92(2 + 9/2) 



(B64) 



(B65) 



(1 + a2 + /2)3 



(B66) 
(B67) 



We provide the ensuing fixed points in table llXi fnote they include the EW and KPZ fixed points), while their stability 
properties are given in table [X] 



Name 


a2 


/2 


92 


z 


a 


/3 


EW 











2 


1/2 


1/4 


KPZ 





0.07583 


0.4642 


1.54 


0.46 


0.299 


FP4,1 


-1.1946 


-1/6 


-0.09436 


1 


1 


1 



TABLE IX: Fixed points of (|B65|l - (fB67l ) in the (a2,/2,ff2) plane. 



Name 


Ai 


A2 


A3 


Stability 


EW 


1 


-2 


1 


Saddle 


KPZ 


0.54 


-2.90 


-1.05 


Saddle 


FP4,1 


88.4 


7.32 


-0.24 


Saddle 



TABLE X: Eigenvalues of the fixed points of table IIXI 



Finally, a further informative set of coupling parameters is 



03 = 



C4A3 



/3 = 



C2A 



53 



A^BqA^ 
4^C3 ' 



(B68) 
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whose flow is 



das „ , 53 



dl 



52ai - 2hal + 4/|a3 + 2/| - ISOa^ - 48a3/3 + 7a3 - 6/| - /a (1 + 03 + /a)" 
.93 (3a3 + h) 



(B69) 
(B70) 

(B71) 

.(l + a3 + /3r 

In this parameter space the fixed points and their exponents are in table IXll and the associated eigenvalues in table 



■ 93 



{l + as + fsY 

93 



Name 


ai 


/i 




z 


a 


/3 


Smooth 











1 








Galilean 








2 


1 


1 


1 


FP4,1 


0.1395 


-0.8371 


0.05534 


1 


1 


1 



TABLE XI: Fixed points of (fB69l) - (lB7T|) in the (03, /a, 53) plane. 



IXIII The Fp4^i fixed point is obtained by solving equation 374a;^ — 2a; — 7 = and using the only admissible solution 
xi = (1 + \/2^)/374 for the components of vector {xi, -6x1, L{xi)), where L{x) = (167238- 1184895a;)/34969. 



Name 


Ai 


A2 


A3 


Stability 


Smooth 
Galilean 

Fp4,l 


-1 
2 
1.09 


-3 

4^*4 

2.28 


-2 

3 - 

4^*4 

-1.98 


Stable 
Unstable 
Saddle 



TABLE XII: Eigenvalues of the fixed points of table IXj 

As done for the |fcp term, we can show the irrelevance of the fc"* "surface diffusion" term, compared to the 
"surface tension" one for C2, C4 > 0. After rescaling of equation the ratio between these two terms is equal to 
C2h'^ /Ci. In the presence of surface diffusion, we obtain several fixed points which are not present in its absence. 
However, the only fixed point with exponents different from those found in section ITVl is Fp4.o, i.e. the linear MBE 
fixed point, while Fp4_i has the same exponents of the Galilean fixed point, z = 1 and a — 1. The stability matrices 
show the latter fixed point has one stable direction, while the Galilean and Fp4^o fixed points are both unstable. The 
KPZ fixed point has two stable directions and becomes the only stable fixed point of the DRG flow if we restrict 
our analysis to the Ci = case, i.e. 04 = 0. Finally, as for the \k\^ term, surface diffusion does not introduce any 
unexpected behavior in case Ci ^ and C2 > 0, hence we conclude that the C4 term is required only when C2 < 0, 
in order to to have a dispersion relation that is well posed in the fc — >■ 00 limit. 

In summary, within the limitations of the DRG technique, we have shown that the main conclusions obtained for a 
linear dispersion relation of the k — k"^ type also apply in the presence of higher order polynomial contributions, as long 
as they have the proper stabilizing signs. Only when the fc^ surface tension term is destabilizing, such higher order 
terms are required for physical consistency. A similar conclusion is expected for /i 7^ 1 values. This irrelevance (and 
also the self-consistency of the DRG calculations) qualifies equation ([55)1 as a faithful representative of the scaling 
behavior of as employed in the main text. 
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FIG. 1: Position of the Galilean fixed point on the (/, g) plane for d = 1 (left panel) and d — 2 (right panel). Arrows and colors 
indicate, respectively, the displacement of this fixed point for increasing fi, and its stability. The intervals of /i specified near 
the solid lines provide the values of ^ at which stability changes for this fixed point. 
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FIG. 2: Real parts of the eigenvalues of the linear stability matrix of the Galilean fixed point for d = 1, as functions of fi. The 
inset is a zoom of the region 4/3 < fi < 5/3 within which the fixed point changes stability from stable to saddle a,t fi — 3/2. 
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FIG. 3: Numerical integration of equations (|41|l - (|42| l for d — 1 and different values of tfie exponent /x: (a) ^ — 1/4, (b) jj, — 1/2, 
(c) i-i — I, (d) n — 1.15, (e) n = 1.4, (f) jj, — 3/2, (g) /i — 1.6, (h) jj, = 1.8. Tlie values of the coupling variables at the different 
fixed points (in the graphs, E stands for Edwards- Wilkinson, K for KPZ, S for Smooth, and G for Galilean) and the associated 
critical exponents are reported in table H] The color of each fixed point represent its stability (red stands for unstable, blue for 
saddle, and green for stable). Note that the Galilean fixed point changes with /i, merging with the KPZ fixed point at jj. — 3/2. 
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FIG. 4: Stability of the fixed points for d = 1,2 and d > 2. As before the red color stands for unstable, blue for saddle and 
green for stable fixed point. 
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FIG. 5: Real parts of the eigenvalues of the linear stability matrix ^(G) for the Galilean fixed point for d = 2, as functions of 
/i. Note the different scales used in the left [for Re (v\j\ and the right [for Re (v2)] vertical axes of the graph. 
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(a) ^ (6) 




FIG. 6: Numerical integration of equations (|41|l - (|42p for d = 2 and different values of the exponent /x: (a) fj, — 1/2, (b) fi = 1, 
(c) fi — 1.3, (d) fj. — 3/2. The values of the coupling variables at the different fixed points (in the graphs, E stands for 
Edwards- Wilkinson, K for KPZ, S for Smooth and G for Galilean) and the associated critical exponents are reported in table 
im The color of each fixed point represent its stability (red stand for unstable, blue for saddle and green for stable). Note that 
the Galilean fixed point is not present for fi > 3/4. 
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FIG. 7: Real parts of the eigenvalues of the linear stability matrix S{G) for the Galilean fixed point for d = 3, as functions of 
/i. Note the different scales used in the left (for v\) and the right (for «2) vertical axes of the graph. 




FIG. 8: Numerical integration of equations (j4T}-(j42j for d = 3 and different values of the exponent [i: (a) [i — 1/2, (b) jj, — 3/2. 
The values of the coupling variables at the different fixed points (in the graphs E stands for Edwards- Wilkinson, R for the fixed 
point related to the roughening transition, S for Smooth and G for Galilean) and the associated critical exponents are reported 
in table Hill The color of each fixed point represent its stability (red stand for unstable, blue for saddle and green for stable). 
Note that the Galilean fixed point is not present for fj, — 3/2. 




FIG. 9: Qualitative behavior of the DRG flow for different values of fi and the substrate dimension d. The color of each fixed 
point represent its stability (red stand for unstable, blue for saddle and green for stable). The shaded region of the plane 
corresponds to the values of fj, (as function of d) within the Galilean fixed point is not defined. 



